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ABSTRACT 


I 


1 


An  experiment  Is  described  which  attempts  to  rn.ea.sure  the 
Isoplanatic  patch  size  of  the  atmosphere  as  a function  of  vari- 
ous parameters.  Technique  constraints  and  the  data  analysis 
procedures  used  are  described  and  results  are  discussed.  It  is 
shown  that  the  constraints  in  data  collection  impaired  the  data 
analysis  and  prevented  a measurement  of  the  isoplanatic  patch 
size.  Recommendations  are  given  for  getting  better  data  in  the 
future . 
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1.0 


SUMMARY 


The  purpose  of  this  project  is  to  estimate  the  isoplana- 
tic  patch  size  of  the  atmosphere  as  a function  of  elevation  angle# 
meteorological  conditions#  and  possibly  o :her  variai'les.  The 
concept  of  the  isoplanatic  patch  is  discussed  in  reference  1. 
Briefly#  a region  of  the  atmosphere  is  termed  an  isoplanatic 
patch  if  the  optical  characteristics  of  the  atmosphere  are  essen- 
tially uniform  oyer  the  region. 

The  method  of  measurement  is  as  follows*.  Pairs  of  stars 
are  selected  which  have  separations  ranging  from  a few  seconds 
of  arc  to  over  200  seconds  of  arc.  Since  stars  are  approximately 
point  sources#  the  image  of  each  star  is  the  point-spread  function 
(PSP)  of  the  atinosphere/telescope  channel  through  which  it  is 
viewed.  If  the  two  images  are  alike#  they  lie  within  an  isoplana- 
tic region.  The  test  of  similarity  of  images  is  the  degree  of 
success  'With  which  one  image  can  be  used  to  restore  th.e  other 
image  by  deconvolution.  The  quality  of  restoration  is  measured 
by  the  width  of  the  image  before-  and  after  processing.  By  pro- 
cessing a sufficient  number  of  star-pair  images  and  obser\'lng 
the  degree  of  similarity  as  a function  of  separation#  it  should 
be  possible  to  estimate  the  probability  that  a patch  of  a given 
size  will  be  isoplanatic.  The  ix>le  of  Federal  Scientific  Corpora- 
tion in  tills  project  was  to  develop  data-reductlon  procedures  for 
evaluating  the  input  data#  doing  the  image  restoration#  and  evalu- 
ating the  success  of  the  restoration#  and  to  perform  the  reduction. 
In  addition#  we  have  presented  some  basic  relationships  which 
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j give  some  Idea  of  the  limitations  of  the  restoration  process  and 

‘ the  tradeoffs  peimltted.  This  background  helps  in  Judging  the 

« 

i restoration  results  Intelligently. 

I • In  analyzing  the  results  of  the  experiment,  we  determined 

I that  the  data  provided  to  us  for  processing  were  unsatisfactory, 

I chiefly  because  the  star  pairs  had  been  photographed  with  too 

long  an  exposure,  but  possibly  also  because  of  scattering  of 
light  within  the  emulsion  of  the  film.  As  a result,  most  of  the 
Images  tended  toward  the  Gaussian  shape,  regardless  of  the  | 

separation  between  the  stars,  and  it  v;as  Impossible  to  tell  ; 

whether  the  two  stars  had  been  Inside  an  Isoplanatlc  patch  or  not.  | 

Further  research  should  be  done  using  shorter  exposures 
and.  If  possible,  thinner  emulsions.  It  is  possible  that  the 
range  of  available  photographic  media  for  this  purpose  has  not 
been  fully  explored.  Alternate  imaging  methods  (e.g.,  vidlcons) 
should  also  be  Investigated. 


f 

i 

i- 

1 


2.0  THEORETICAL  STUDISS 

2 . 1 V/ldth  Improvement 

In  the  following  discussion,  we  assiome  familiarity  with 
the  general  idea  of  image  restoration  by  division  of  Fourier 
transforms.  (Refs.  2,  3>  ^.) 

Let  the  point-spread  function  (PSF)  of  the  atmosphere/ 
telescope  channel  be  h(x,y),  the  blurred  image  b(x,y)  and  the 
ideal  image  p(x,y);  let  capitals  denote  P’ourier  transforms  and 
F indicate  Inverse  transformation.  Then,  in  the  absence  of 
noise,  the  ideal  image  can  be  recovered  by  means  of  the  relation, 

fiM  = (S.  ?;]  (1) 

assuming  that  at  zeroes  of  H,  the  limit  of  Ey^H  exists.  In  prac- 
tice, both  b and  h are  corrupted  by  noise.  This  can  be  haze  or 
other  s<;atterecl  light  in  the  atmospheres  fogging  in  the  film, 
grain,  or  sampling  noise.  Some  of  tnis  noise  is  additive  and 
some  not.  If  we  assume  it  to  be  additive,  we  can  write 
d(x,y)  - b(x,y)  + n^(x,y) 
k(x,y)  = h(x,y)  + rij^(x,y) 

Because  of  this  noise,  deconvolution  by  transform  methods  gives 
us  only  an  estimate  of  p: 

,-S  (^,y)  “ 

In  practice,  and  are  nearly  uniform  over  the  transform  do- 

\ 

main,  while  D and  H fall  off  rapidly  with  increasing  spatial  1 

frequency.  As  a result,  the  high-frequency  regions  of  D/K  are 


3 


3 ^ 


(2) 
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extremely  noisy  and  make  little  contribution  to  Jo  , It  is  cus- 
tomary to  deal  with  this  problem  by  some  kind  of  low-pass  fil- 
tering of  D/K  (see,  for  example,  reference  3).  It  would  be  nice 
if  we  could  design  this  filter  to  optimize  the  approximation  of 
to  p in  a mean-square  sense.  Analysis  indicates,  however,  that 
such  a filter  has  the  effect  of  cancelling  the  inverse  filter 
and  replacing  it  with  a minimum  MSE  filter.  Since  the  object  of 
the  present  experiment  1s  to  use  success  in  Inverse  filtering  as 
a measure  of  Isoplanaticlty , this  replacement  is  not  desirable. 
Instead,  we  use  a simple  rectangular  filter: 


O oth^rn'/S'^ 


(3) 


Then  the  restored  image  1s 


P (^,y)  = ^ 


-/ 


^ /-A/^ 

Experience  has  shown  that  it  is  nearly  always  possible  to  reduce 
the  apparent  noise  in  ^ to  any  desired  degree  by  choosing  R suf- 
ficiently small.  The  resulting  image,  however,  is  convolved  with 
g^,  the  Inverse  transfomi  of  the  filtering  function,  and  as  R 
decreases,  the  blurring  effect  of  may  become  worse  than  that 
of  h,  which  we  are  trying  to  remove.  Indeed,  if  R is  chosen  so 
that  the  transforms  are  virtually  noise-free  Inside  R,  then  the 
operation 


// 


// 


//  A// 


essentially  exchanges  a new  PSP  (Sq)  for  the  old  one  (h).  Unless 
Sq  Is  Jn  some  sense  "tighter"  than  h,  this  Is  futile.  In  light 
of  thlSj  It  is  of  Interest  to  ask^  for  what  value  of  R is  the 
blurring  effect  of  gQ  worse  than  that  of  h? 

The  answer  to  this  question  depends  on  the  shape  of  h and 
on  the  way  in  which  we  measure  the  blurring  effect,  or  "width," 
of  the  PSP.  For  this  analysis,  we  assume  h Gaussian  and  centered 
at  the  origin,  and  measure  width  by  the  volume/height  ratio. 


= ZS2.<3S. \0l 

^Co.o) 

(This  is  similar  to  the  principal  width  measure  used  in  the 
present  project.) 

If  we  assume  circular  symmetry,  it  is  convenient  to  use 
Hankel  transforms.  Let  r = radial  distance  and  w = radial  fre- 
quency. Then  we  have  the  following  transform  pairs: 


1 
2 I 


Filter 


iV  < 


J,  (kr) 


^ I ry/  - 3 \ - 

2 ^ 

(y  = i3  the  Airy  pattern).  Using  an  N-polnt  discrete 

finite  Fourier  transform,  however,  the  sample  spacing  Is  not  in 
radians  but  in  units  of  l/N  cycles.  If  we  let = radial  fre- 
quency In  samples,  /o  = and  the  transform  pairs  become* 

♦These  relations  are  approximate,  since  they  do  not  take  Into 
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PSF 


where 


i'« 


Filter 


\ 2 Jf  (^2  Jr  a-/aJ  ) 

— ' 7T^ 


y ^ 


The  width  of  the  Gaussian  PSF  Is 


(( 

c>  othiPrK'/se’  '• 


% --F.7J-az^  . 

" " Zfra'^^ 


and  the  width  of  the  Airy  pattern  is 


(10) 


^ 1^/e  ^ 

/ 


(11) 


Then  the  "break-even"  point  is  the  value  of  Ri  for  which  W = W,  : 

g h 


If  noise  conditions  force  us  to  filter  at  some  radius  , then 
the  expected  Improvement  In  image  spread  Is 


(12) 


(13) 


■ 

i. 


2.2  Noise  Considerations 

In  practice,  one  tries  various  values  of  the  filter 
radius,  Rj,  until  the  best  compromise  Is  obtained  between  en- 
hancement and  noise  level.  The  limits  that  govern  this  compromise 


account  the  wrap-around  and  aliasing  associated  with  finite 
sampled  domains.  With  suitable  choice  of  Image  size  and  sampling 
rate,  however,  the  approximations  are  close  enough. 
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can  be  estimated  in  the  following  way:  The  blurring  of  an  image 

can  be  considered  the  result  of  passing  it  through  a channel 
whose  OTF,  ) wipes  out  high-spatial-frequency  components 

of  the  image.  From  elementary  information  theory,  the  maximum 
information  the  blurred  image  can  contain  is 

h = ^ Gits)  (14) 

v,’here 


OO  CO 

LL\ 

(2 

I ^ (o,c) 

2 

and 

(16) 

is  the  signal-to-nolso  ratio  in  the  image  and  and  y^  are 
the  dimensions  of  the  image. 

Now  there  are  two  Inputs  to  the  restoration  process  — 
the  blurred  image  and  the  PSF  --  and  the  information  in  the  re- 
stored image  can  come  only  from  these  inputs  If  the  information 
in  the  PSF  is  = ®h^o^o^h  information  in  the  restored 

image  is  Ip  = BpXgyQkp  (where  Bp,  and  kp  are  defined 
analogously  to  equations  15  and  l6) , then 

(IT) 

With  equality  only  if  .:he  restoration  process  takes  advantage  of 
all  the  available  information.  If  we  have  Isoplanatlclty,  so 
that  b(x,  y)  and  h(x,  y)  go  through  identical  channels,  then  B^  = 
B^j  hence  we  get  the  following  tradeoff  between  enhancement  and 


noise : 


< 


(18) 


(/  -f  S^)  + (/  y-  S'j,) 

/oa.  (/  y-  sy 


In  the  case  of  star  pairs.  If  the  stars  are  of  equal  magnitude 
so  that  = S^,  and  if  we  assume  that  the  narrowing  of  the  image 

is  roughly  proportional  to  the  widening  of  the  bandwidth,  then 
equation  l8  says  that  we  can  get,  at  best,  a 2-to-l  Improvement 
in  image  size  if  we  will  accept  no  degradation  of  signal-to- 
noise  ratio. 

If  we  return  to  our  Gaussian  assumption,  we  can  be  more 
specific.  For  the  Gaussian  PSF,  B,  foP  the  filter  of 

rt  // 

equation  9,  Bf  . Then,  using  the  definition  of  equation  13, 


< 


^ C ^ ^ f / s^') 

1 (/-^  ~^)  (^9) 

These  results  are  generally  consistent  with  the  experi- 
mental results.  In  particular,  for  stars  with  approximately 
Gaussian  Images,  the  actual  restoration  results  closely  parallel 
equation  ±9,  in  spite  of  the  rather  crude  method  used  for  esti- 
mating slgnal-to-nolse  ratio.  Among  the  representative  results 
singled  out  in  Section  k,  the  only  stars  whose  results  exceed 
equation  19  are  those  whose  Images  are  clearly  non-Gausslan.  For 
the  other  stars,  the  results  are  almost  uniformly  about  80^  of 
the  maximum  predicted  by  equation  19. 
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3 . 0 EXFSRIMllNTAL  METHOD 

3 . 1 Data  Acquisition 

Data  were  obtained  by  photographing  star  pairs  at  the 
f/10  focus  of  the  ^8-lnch  telescope  at  the  Range  Measurements 
Laboratory  (RML),  Patrick  Air  Force  Base.  A list  of  the  star 
pairs  obser-ved  is  given  in  Table  I.  The  magnitudes  range  from 

2.1  to  6 and  the  separations  fr-oin  4.7  to  28l  seconds  of  arc. 

The  exposure  time  In  most  cases  -was  1/'125  of  a second.  No  effort 
was  made  to  photograph  star  pairs  under  "good  seeing"  conditions. 
Approximately  six  photographs,  spaced  about  one  second  apart, 
were  taken  of  each  pair  at  several  elevation  angles.  The  photo- 
graphs were  made  using  Kodak  Tri-X  film,  which  v;as  developed  for 
an  ASA  of  600  to  800.  One  frame  of  each  batch  of  film  v;as  ex- 
posed to  a wedge  filter;  this  image  was  used  to  remove  the  H & D 
curve  from  the  data.  The  final  result  thus  indicated  intensity 
instead  of  density. 

For  each  star  pair,  four  to  six  best-looking  images  were 
selected  for  digitization.  The  star  images  were  digitized  from 
the  negative  by  means  of  the  laser  iraage-processlne  scanner  at 
RML,  using  a 5"mlcron  spot  size  (corresponding  to  a resolution 
of  0.85  second  of  arc)  and  an  8-bit  word  length.  The  digitiza- 
tion generated  a 128-x-128-point  matrix.  This  was  then  edited 
into  two  64-x-64-polnt  matrices,  each  containing  one  star  image. 
Matrices  wr>re  written  on  magnetic  tape  and  sent  to  FSG  for  pro- 
cessing. A single  tape  typically  contained  several  files,  each 
one  consisting  of  three  records:  a descriptive  header  record  and 
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records  containing  the  edited  images  of  two  stars.  Because  errors 
sometimes  occurred  in  editing-out  individual  star  Images # files 
containing  the  original  128-x-128-point  matrices  were  also  in- 
cluded on  the  tape. 

The  header  record  Included 

- File  number 

- Sequence  number 

- Matrix  size 

- Star  identification,  right  ascension, 
declination,  separatxon,  and  visual 

magnitude  (of  brighter  star)  |j 

- Azimuth  and  elevation  at  time  of  exposure 

- Time  of  exposu3?e 

- Date 

- Duration  of  exposure 

- Coordinates  of  matrix  origin 

Each  image  was  written  as  a set  of  records,  one  record  for  each 
row  of  the  image  matrix. 

3.2  Lata-Frocesslng  Procedures 

Each  star  tape,  upon  arrival,  was  given  a number,  and  a 
table  of  contents  was  generated  by  listing  its  headers.  From  the 
table  of  contents,  the  files  containing  64-x-64-polnt  arrays  were 
identified.  These  files  were  then  processed  by  two  programs:  the 
preliminary  program  and  the  image-processing  program.  These  will 
now  be  described. 

The  preliminary  program  reads  the  files,  plots  the  images 
contained,  and  Fourler-transforms  the  images.  In  reading  the 
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images,  the  program  converts  the  f3xed-point  8-blt  samples  into 
32-bit  floating-point  numbers,  scaling  them  so  that  the  sample 
range  of  0 to  255  converted  into  data  values  ranging  from  0 
to  2.55,  The  image  matrix  used  in  processing  is  a 64-by-64 
array  of  complex  numbers;  the  program  takes  the  floated  sanple 
values  and  inserts  them  in  the  real  pa.rt  of  this  complex  array. 
Thus,  for  example,  a sample  value  of  "(I'  (hexadecimal)  ends  up  as 
the  complex  number  (l.27>  0.).  The  Images  are  plotted;  these 
plots  enable  us  to  Judge  image  quality  and  to  reject  any  im.ages 
containing  obvious  errors.  For  each  pair  of  images,  the  program 
writes  a set  of  output  data  containing  the  images  and  their 
transforms.  The  sequence  is: 

Header  (croied  from  input  tape) 

First  ima<i,.i 

Transform  of  first  image 

Header 

Second  image 

Transform  of  second  image 

The  program  provides  means  for  specifying  the  order  in  which  the 
images  are  processed.  In  most  cases,  the  smaller  image  of  the 
pair  was  selected  first.  (This  is  because  th^'  first  image  is 
used  as  the  reference  PSP  in  the  image-processing  program,  and  it 
is  not  reasonable  to  expect  an  image  to  be  smaller  than  the  PSF 
that  gave  rise  to  it.) 

The  image-processing  program  does  the  remainder  of  the 
data  reduction.  Thjs  consists  of  the  following  steps: 


i 

I 

j 


li, 

\ 

! 


I 

I 
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Evaluate  first  Image 
Evaluate  second  image 
Filter  second  image 
De-convolve 

Evaluate  resultant  image 

The  image-evaluation,  filtering,  and  de-convolution  routines  will 
now  be  described. 

Image  Evaluation.  These  routines  compute  image  centroid, 
maximum  amplitude,  slgnal-to-noise  ratios,  and  four  measures  of 
image  width.  The  maximum  and  centroid  computations  are  trivial 
and  will  not  be  described  here.  Signal-to-nolse  ratio  is  esti- 
mated by  taking  the  maximiun  amplitude  as  the  signal  and  the  stan- 
dard deviation  of  the  values  about  the  periphery  of  the  matrix  as 
the  noise.  This  is  an  approximate  estimate;  it  rests  on  the  as- 
sumption that  thf  Lm.age  is  narrow  and  approximately  centered  in 
the  matrix,  so  that  the  values  at  the  edges  are  m.ostly  noise 
These  assumptions  are  usually  justified,  but  occasionally,  v;hen 
the  linage  is  unusually  big  or  is  badly  off-center,  the  noise 
estimate  goes  bad  and  the  measure  is  worthless.  It  Is  usually 
easy  to  detect  these  occurrences,  however. 

The  image-width  measures  are: 
a . Squal-voliurne  area 


/V.  = 


Oj')) 


b.  Second-moment  area  ^ 

_ __  ^ 


= — 


5 


12 


c. 


^ O,-^)  — -707  MiZ-K 


Measures  (£)  through  (c)  are  computed  by  a single  subi*outine.  In 
the  case  of  a restored  image,  the  intensity  can  sometimes  go 
negative  (reference  4);  this  creates  problems  for  the  first 
three  measures  and  maRes  their  interpretation  difficult.  To 
avoid  this,  the  program  sums  only  over  the  central  lobe  of  the 
image,  stopping  at  the  contour  at  which  the  amplitude  goes  to 
zero.  This  contour  is  found  by  a sea‘»^ch  program  which  tags  all 
pixels  belonging  to  the  central  lobe.  The  result  of  this  re- 
striction is  that  the  width  estimates  for  the  restored  image  tend 
to  err  by  being  too  optimistic.  That  is,  the  numerical  measure 
of  image  narrowing  seems  high  compared  to  what  we  would  expect 
from  visually  ccmpai’ing  the  before-and-after  plots.  Without  this 
m.aln-lobe-  restriction,  however,  the  width  estimates  for  the  re- 
stored images  err  grossly  in  the  pessimistic  direction,  and  also 
show  little  or  no  coxislstency . V/e  decided  that  the  main-lobe 
measures  were  superior,  even  with  the  bias. 

Image  Filtering.  The  filter  function  zeroes  cut  all 
values  in  the  transfom  for  which  the  noise  is  dominant.  We 
assume  circular  syimnetry,  since  this  is  the  best  assumption  in 
the  absence  of  ^ priori  knowledge  about  the  image  spectrum.  The 

13 
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filtering  function  used  shown  in  Fig.  1. 

/,  A < ^,/£ 

f'  r TT  \ ^ y 

?"y ^ 3^,/2 

o,  r>  3 e,/Z 


\ e,  3e, 

Pig.  i.  Spectrum  low-pass  filtering  function 

R]_  Is  referred  to  as  the  truncation  radius.  In  choosing 
Rj^  there  is  a trade-off  betv.een  image-v;ldth  reduction  and  slgnal- 
to-nolse  ratio,  as  discussed  above  (Section  2).  The  cosine- 
squared  rolloff  Is  used  to  suppress  sidelobes  In  the  restored 
image. 

The  filtering  program  computes  t^he  filtering  function 
for  a specified  value  of  R-,  and  multiplies  it  by  the  transform 
of  the  image. 

De-convolution.  Conceptually,  the  de-convolutlon  pregrarn 

t 

is  trivial:  it  simply  divides  one  Fourier  transform  by  another. 

(The  first  image  is  taken  as  the  reference  PSF  and  its  transform 
(already  on  the  working  tape)  as  the  denominator  transform;  the 
second  image  is  the  one  which  we  attempt  to  restore  and  its 
transform  is  the  numerator.)  In  practice,  the  program  must  be 
able  to  detect  unreasonable  results  and  try  to  replace  them  with 
reasonable.*  ones.  The  unreasonable  results  occur  at  low  S/N  levels, 
where  noise  drives  the  denominato.r  transform  toward  zero.  The 
program  attempts  to  r'*medy  this  state  of  affairs  by  interpolating 
from  neighboring  points,  assuming  that  these  are  not  as  badly 


/.o 


/(>;  = 
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affected  by  noise. 


V/hat  happens  is  this:  Typically,  the  S/N  ratio  near  the 

origin  o^^  the  Fourier  domain  is  high  and  there  is  no  trouble.  As 
we  move  away  from  the  origin,  Ihe  signal,  in  both  the  numierator 
and  the  denominatoi',  drops  rapidl5».  Presently,  we  encounter  iso- 
lated points  v;here  the  noise  has  driven  the  denominator  close  to 
zero.  In  this  region,  the  program  attempts  to  patch  up  the  de- 
nominator by  interpolation.  As  we  move  into  lower  and  lower  S/N 
regions,  the  number  of  such  points  increases  until  their  density 
is  so  great  that  interpolation  does  no  good.  This  "hopeless" 
region  is  supposed  to  be  eliminated  by  the  image-filtering  pro- 
cess describe  •,  previously.  (When  the  numerator  is  zero,  the 
progi'am  looks  no  further  but  immediately  sets  the  quotient  to 
zero  ap  that  point.)  Thus,  the  purpose  of  interpolation  is  to 
Improve  results  in  the  borderline  region  where  only  isolated 
points  are  seriously  corrupted  by  noise. 

If  two  stars  are  viewed  through  a perfectly  Isoplanatlc 
system,  their  PSFs  will  differ  only  by  a constant  multiplier, 
corresponding  to  their  difference  in  magnitude.  Their  Fourier 
transforms  will  similarly  differ,  and  the  deconvolutlo  i quotient 
will  be  a constant  equal  to  the  ratio  of  the  two  transforms. 

Hence  tills  ratio  represents  the  ideal  result,  and  "unreasonable" 
values  are  values  which  differ  grossly  from  this  ideal. 

In  practice,  the  ratio  is  not  constant.  The  program,  ap- 
proximates the  ratio  by  dividing  the  Integral  of  the  image  by 
that  of  the  reference  PSF.  (These  numbers  are  already  available 
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since  they  are  the  origin  values  of  the  respective  transforms.) 
An  "unreasonable"  quotient  point  is  defined  as  one  whose  magni- 
tude is  greater  than  two  times  this  ratio.  (This  threshold  Is 
arbitrarily  chosen.) 

When  an  "unreasonable"  point  is  detected j the  program 
does  a four-point  interpolation  of  the  denominator.  That  is,, 
if  Che  initial  denominator  point  is  d(l#j).  then  the  program 
computes 

d 'Oj 0^4  +■  d 0, J-*- 0 ] 

and  tr-ies  again.  If  the  result  is  still  unreasonable^  the  pro- 
gram throws  In  four  more  points  — 


^ 0, 0 ~ j:  Od)  -'“k  \_^ d O'-/,  y-^/j 

/■d {(  i /^  j -/)  -f-  d (i-f  /,  yd-/)'\ 

and  tries  a third  time.  If  the  quotient  is  still  unacceptable;, 
then  its  magnitude  is  forced  to  the  ideal  value  and  the  point  1 
flagged  for  further  attention. 

After  the  process  is  complete,  the  flagged  values  are 
examined.  A four-point  Interpolation.,  similar  to  that  describe 
above  for  the  denominator,  is  attempted.,  but  only  if  the  four 
reference  points  are  uaflagged.  (This  technique  s not  as  suc- 
cessful as  might  be  expected,  bev-'ause  unsatisfactory  points  cen 
to  cluster. ) 


^.0 


EXPERIMENTAL  RESULTS 


The  results  obtained  for  all  star-pairs  processed  are 
presented  in  Tables  2 through  10  at  the  end  of  this  section. 
Following  these  tables ^ plots  for  selected  pairs  are  shown  In 
Figures  2 through  20.  These  plots  are  also  referenced  at  the 
corresponding  locations  in  the  tables.  Most  of  the  table  entries 
are  self-explanatory.  "Sequence  No."  refers  to  the  tape  and  file 
numbers  which  Identify  the  Image.  In  most  cases,  the  smaller 
image  was  used  as  the  PSF,  but  occasionally  the  Images  were  in- 
terchanged, particularly  if  they  v;ere  approximately  the  same 
size,  to  see  whether  better  results  were  obtained  that  way.  "Im- 
provement" is  given  as  the  ratio  of  v;idth  before  processing  to 
width  after  processing:  thus,  the  larger  the  ratio  the  better 
the  Improvement. 

The  descriptions  under  "Appearance"  are  subjective  evalua- 
tions of  restored-image  quality  as  Judged  irom  looking  at  the  3~D 
plots.  They  are  Intended  to  supplement,  and  in  some  cases  to 
correct,  the  numerical  entries.  As  we  mentioned  in  the  preceding 
section,  the  numerical  results  of  linage  evaluation  are  not  always 
consistent.  Although  descriptions  such  as  "very  good"  and  "fair" 
are  necessaidly  vapjue,  they  are  at  least  not  misleading.  In  ad- 
dition, their  slgnlfic-ance  is  clarified  by  the  3-D  plots  supplied 
for  selected  stars. 

Restoration  with  any  given  star  pair  v/ac  largely  a cut- 
and-try  procedure  using  different  values  of  the  filter  I’adlus, 

Rl.  Where  practical,  we  tried  to  fjnd  a value  giving  the  best 
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trade-off  between  Image  size  and  noise  level.  Because  available 
computer  time  was  limited,  we  were  not  always  able  to  try  as  many 
values  as  might  be  desired.  In  some  cases,  we  quickly  found  that 
no  satisfactory  trade-off  was  to  be  found  and  the  search  was 


abandoned. 


Tables  9 and  10  present  a selected  subset  of  the  processing 
results.  (The  rule  for  selection  Is  to  choose  the  result  for 
which  the  estimated  slgnal-to-noise  ratio  In  the  restored  Image 
is  closest  to  30  dE. ) Inspection  of  these  tables  shows  the 
following: 

1.  Width  enhancement,  using  the  equivalent-area  measure 
(Wl)  averages  I.71.  (A  ratio  of  1.00  indicates  no  improvement 
at  all. ) 

2.  There  is  no  correlation  between  enhancement  and 
separation. 

3.  There  Is  no  correlation  between  enhancement  and 
elevation  angle. 

lixamlnatlon  of  the  star  data  reveals  two  significant 
facts.  First,  nearly  all  star  Images  look  Gaussian.  This  Is 
Inconsistent  with  v;hat  we  had  been  led  to  expect:  Images  degraded 
by  turbulence  were  supposed  to  be  "broken  up"  and  "torn  apart." 
Im*age  tear  was  never  seen,  however.  The  Images  In  Figures  2 and 
3 are  representative  of  what  was  obtained  most  of  the  time.  In 
fact,  most  of  the  Images  of  most  of  the  stars  looked  alike.  To 
test  this,  we  tried  restoring  two  unrelated  Images.  This  ex- 
periment was  run  twice,  first  with  two  Images  from  star-pair  48 
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taken  some  seconds  apai't;  second  with  two  Images  taken  two  months 
apart.  The  results j shown  in  Table  11  and  Figures  l8  and  19,  are 
roughly  comparable  to  the  average  enhancements  obtained  in  the 
main  experiment.  Finally,  enhancement  was  attempted  with  a com- 
puter-generated Gaussian  pulse,  tailored  to  fit  the  star  image. 

The  result  (Fig.  20  and  Table  11)  was  slightly  tetter  than  the 
average,  and  distinctly  better  than  what  we  had  achieved  using 
the  star's  own  companion  image  (Fig.  5). 

If  the  images  are  Gaussian,  then  it  should  be  possible  to 
check  the  results  in  Table  9 against  equation  19.  In  spite  of 
the  rather  crude  method  used  for  estimating  the  slgnal-to-noise 
ratio,  the  results  are  generally  consistent.  The  exceptions  are 
stars  whose  images  are  visibly  non-Gaussian,  notably  the  1-msec 
exposure  of  Star  ^}0. 

The  other  unexpected  result  is  that  the  images  are  not  of 
equal  size.  The  brighter  star  always  has  a wider  image,  and  the 
difference  in  width  is  roughly  proportional  to  the  difference  in 
magnitude.  In  a linear,  isoplanatic  system,  the  star  images  can 
be  expected  to  differ  by  noise  and  by  a multiplicative  constant. 
Neither  of  these  differences  can  cause  the  differences  in  image 
size  seen  here. 

One  possible  explanation  for  the  Gaussian  shape  is  V lurrlng 
resulting  from  too  long  an  exposure.  (See,  for  example,  referen- 
ces 3 and  ) This  possibility  is  suggested  by  the  extraordi- 
narily good  results  obtained  with  the  one  1-msec  exposure  in- 
vestigated. Star  ^10,  sequence  #7/7  and  7/8,  shown  in  Table  9 
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and  in  Figs.  6,  7,  8,  and  9,  yielded  the  best  overall  enhance- 
ment of  all  the  Images  processed.  The  evidence  of  a single  image 
cannot  be  conclusive,  of  course,  but  it  is  suggestive,  especially 
since  the  Images  come  closer  to  showing  Image  tear  than  any  others. 

The  unequal  size  is  most  likely  the  result  of  non-lineari- 
ties in  recording  the  received  Images.  Either  the  exposure  was 
so  weak  bhat  the  images  fell  at  the  non-linear  "toe"  of  the  H & 

D curve  for  Trl-X,  or  else  the  logarithmic  H & D curve  itself 
was  Improperly  removed.  Another  possibility  is  that  ll^t  scat- 
ter occurred  within  the  relatively  thick  emulsion  of  the  Trl-X 
film.  (Scatter  would  also  tend  to  make  torn  images  look  Gaussian.) 
This  would  suggest  that  the  exposure  was  too  strong,  however, 
rather  than  too  weak.  Calculations  based  on  visual  magnitude, 
aperture  size,  etc.,  have  not  been  conclusive,  but  the  images, 
as  seen  on  the  film,  do  not  look  dense  enough  for  scatter  to  be 
likely. 
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TABLE  1 STAR-PAIR  DATA  RECEIVED  FROM  RML 


These  data  not  supplied 


TABLE  2.  SWJ^ARY  OF  RESULTS  FOR  STAR-PAIR  #2^ 


Notes 


TABLE  3.  SUMMARY  OF  RESULTS  FOR  STAR-PAIR 


TAELS  4.  3UMHARY  OF  RESULTS  FOR  STAR-PAIR  #32 
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Noteei  1.  Az,  el  eJata  not  supplied  for  these  images 

2.  Image  filtered  only,  not  restored 

3.  These  images  are  illustrated  in  Pigs.  1^, 


TABLE  10  REPRESENTATIVE  RESULTS  IN  ORDER  OF  ELEVATION  ANGLE 
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CONCLUSIONS  AND  RECOMMENDATIONS 

Most  of  the  star-image  pairs  supplied  for  this  study 


showed  no  image  tear  and  tended  to  look  Gaussian.  The  only  ex- 
ception was  the  star-pair  photographed  at  a 1-msec  exposure. 

Images  of  stars  of  unequal  magnitude  were  of  unequal  width. 

These  facts  suggest  that  the  star-image  data  supplied  for  this 
study  were  unsatisfactory. 

Efforts  for  the  immediate  future  should  concentrate  on 
remedying  these  defects.  The  Gaussian  image  appears  to  be  the 
result  of  too  long  an  exposure.  The  unequal  image  size  may  be 
the  result  either  of  scattering  within  the  e'  vision  or  of  under- 
exposed film.  If  it  is  the  latter,  it  will  be  difficult  to  solve 
the  problem  without  going  to  a longer  exposure  and  thus  aggravat- 
ing the  Gausslan-lmage  pi’oblem.  It  may  be  necessary  to  investi- 
gate other  types  of  emulsion  or  other  means  of  recording  the 
image.  Alternatively,  future  experiments  may  have  to  be  restricted 
to  photographing  only  the  brightest  stars. 

As  collateral  studies,  experiments  might  be  run  using 
artificial  sources,  e.g.,  high-intensity  lights  or  lasers,  either 
flown  or  at  ground  level,  especially  if  exposure  constraints  make 
it  difficult  to  find  usable  star-j)airs  covering  a wide  range  of 
separations. 
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